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A method is proposed for calculating the absorption spectrum of a long one- 
dimensional closed-into-a-ring chain with Frenkel's exciton under diagonal disorder. 
This disorder is represented by the hyperbolic singularities of atomic fission. These 
defects are shown to lead to a wing in the exciton zone of a chain without defects. The 
form of the wing does not depend on the relative positions or number of defects and 
its value is proportional to the sum of the amplitudes of the defects. The proposed 
method uses only the continual approximation. 

1. Introduction 

Progress in the theory of translation-symmetric systems is largely due to the possibility of 
using the Bloch function tools. Switching to the Bloch representation makes it possible to 
appreciably simplify and, sometimes, even to solve a variety of problems in the physics of 
translation-symmetric systems. The absence of a similar universal approach in the theory of 
systems without translation symmetry causes many mathematical difficulties in the analysis 
of even the simplest systems of this kind. One such system is a model with the Frenkel 
exciton-type excitations, i.e. a set of two level atoms bound by an interaction capable of 
earring the excitation from one atom to another. One usually explores only the one-exciton 
region of energy spectrum, which corresponds to the presence of a single excited atom. If 
the atoms are arranged symmetrically and have different level splits, the system is called 
diagonally disordered. Two problems of this kind - the Dyson problemB and the Lloyd 
problem! - are two of the very few that allow an exact solution. Approximate methods 
(such as the coherent potential approximation!) were developed for these systems. The 
quality criterion for these methods is often a comparison with data from computer experi- 
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ments. However, it is not usually clear why the approximate calculations agrees (if it does) 
with the numerical calculation. Therefore, in the author's opinion, even abstract models 
of translation-nonsymmetric systems are of interest if a mathematically convincing solution 
scheme can be proposed for them. In the present paper, such an abstract model is proposed. 

We consider the problem of the Frenkel exciton in an infinitely long chain of two- level 
atoms that has been closed into a ring. The interatomic interaction and the atomic fission 
energy as a function of the atom's coordinate are 

w{z) = VeM-\^\), (1) 
and 

^ a 

respectively, i.e. diagonal disorder is introdused in the form of a random number N 
of hyperbolic singularities (hereinafter, hyperbolic defects) which are randomly placed at 
points Rr and have positive amplitudes ar (of constatnt sign). In present paper, we propose 
a method of calculating the absorption spectrum of this chain at zero temperature. 

The absorption spectrum of a defectless chain is known to be a singlet. Therfore, changes 
in the spectrum due to defects are easier to see than changes in the density of states. 
It is shown herein that hyperbolic defects give rise to a wing in the exciton zone of a 
defectless chain and that its shape does not depend on the number or relative positions of 

N 

the defects. Also its size is proportional to the sum of the amplitudes of the defects, Yl,cbr- 

r 

This result appears to be somewhat unexpected. The calculation was carried out by the 
method proposed in@, which uses only the continuation approximation (i.e. the replacement 
of the lattice sums by integrals). The efficiency of this approach was substantiated in@ 



2. Scheme of the method 



In oder to calculate the absorption spestrum of the chain discribed in the Introduction, we 
use the method proposed in@. Here we give the results from that paper. If each spectral line 
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is thought of as a Lorentz hne with width S, then the absorption spectrum A{Q) of a chain 
with parameters (1),(2) and end coordinates L and —L can be calculated by the formula 

-L ^ ' 

where ^{z) is defined by 

Here R is the radius of interaction {\\W = 2pVR is the width of the exciton zone of 
the defectless chain, p is the density of the atoms in the chain, e is the atomic fission energy, 
Vt is the energy of the incident light quantum, and E = VL + i5,5>Q. For a chain that 
closed into a ring, the function '^{z) must satisfy the cyclicity conditions 

^(L) = ^(-L) (5) 

d-^{L)/dz = d'^{-L)/dz (6) 

Integrating (4) over z from L to —L and taking into account (5) and (6), we find that 
the integral in (3) can be represented as follows: 

--/^-/«^'-^)- m 

-L ^ ' -L 

L 

Thus the task is redused to calculating the integral / ^{z)dz. The idea is to calculate 

-L 

this integral not along the real axis but over semicircle of radius L in the upper half-plane of 
the complex variable z. This is possible if '^{z) does not have singularities (poles or branch 
points) in the upper half-plane. Next, if e{z) decreases in the above semicircle, we have 
the following relationship: the lager the semicircle radius (chain lenth), the more accurately 
Eq(4) can be solved by perturbation theory on semicircle. As L — cxd, the error tends to 
zero. It is known! that the solution to (4) has no singularities in the upper half-plane of 
complex z if the function is a single-valued and has no poles there. 



Let us show that this is the case for e{z) of form (2), where E has a small positive 
imaginary part. If E is real, the equation E — e{z) = has N real roots Zk{E), k = 1, ...,N . 
If E has small imaginary part term i5, 5 > 0, the roots obtain an increment 



5 — 
^ de{zk)/dz 

Since €{z) defined by (2) is decreasing, all roots are shifted to the lower half-plane and, 
thus, in our case, (^e^{z)) ^^^^ ^^t have poles in the upper half-plane. In the next section, 
we calculate the spectrum by discribed scheme. 

3. Calculation of the spectrum 

Let us represent (4) in the form 

+ g^* + A(z)^ = -{l/ny (8) 
where 



RJ \E 

Here q is the wave vector of the exciton with energy E in the defectless chain, the function 
A(z) describes the defects and decreases in the semicircle of radius L (below we call it large 
semicircle) in the upper half-plane of complex z. Let us find the solution of (8) in the form 
^{z) — ^o{z) + + where the lower index indicates the order of correction A.{z). 

These corrections obey the following equations: 

^'o' + q'^o = -{l/Rf (11) 

+ = -A*o (12) 
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It is not difficult to solve Eqs.(ll) and (12). After that, in the approximation under 
consideration, ^{z) becomes 

L 



+C+ [ e''^' - J —^e'''^ sm[q{z - 0]d^j + 



+C_ I e-'"" ' ^^^^ 



All integrations begin at the right end of the chain and go over the large semicircle. 
Since l^{z) is small there, (as L — > oo), we expect that solution (13) in the large semicircle 
virtually coincides with the exact solution. To find the constants C+ and C_, we use the 
solution (13) and conditions (5) and (6), which is possible because the points ±L lie in the 
large semicircle where this solution is vahd. The result is the following set of equations for 
C+ and C_: 

C+{lM.ct) - 2qsm{qL)} + C^{h{-q,q) - 2qsm{qL)} = h{0,q)/{Rqf (14) 

C+{tIs{q,q) - 2qsm{qL)} + C_{tls{-q,q) + 2qsm{qL)} = il , q) / {Rq)\ (15) 
where 

-L 

I,(a, (3)= I A(Oe*"« cos[/3(L + Ol^C (16) 

L 

-L 

I,{a, /3) = - / A(Oe^"^ sin[/3(L + OH (17) 
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Recall that all integrations are carried out over the large semicircle. We consider q to 
be real and positive. This means that the energy E falls within the exciton zone 
of defectless lattice. Now we calculate integrals like (16), (17) that appear in (14), (15). 
Wherever possible, we proceed to the limit L — > oo. Let us explain, for example, how 



— 1^ 



is calculated. In the large semicircle, only the increasing exponent exp{—iq{L+^)) should 
be retained in sine function. On the other hand, this exponent tends to zero in the lower 
large semicircle. Therefore, 

"/A(OexpHg(L + 0]rfe 

L 

can be calculated over the entire large circle, which is quite easy to do with the help of 
residues. We give the results of calculating the integrals appearing in (14), (15) below: 

ls{0,q) = 7re-^'^^^Res[A(0e-^''«] (18) 

1 

I,{q,q)^-e-^<i^ J A{OdC (19) 

L 

/e(0, q) = TTie"*"^ Yl Res[A(Oe-*''«] (20) 

1 7"- 

Ic{q,q)^^e-^''' j (21) 

L 

Since /c(0, Q') = ils{0,q) and ils{q,q) — Ic{q,q), then C_ = and we do not need the 
remainig integrals from (14), (15). For C+, we obtain the following formula: 

= U — ^ '''^ 

g-igL J A(0ciC-4gsin(gL) 

L 

-L 

For the final calculation of / (^(^) — 1)0?^;, we also need the following integrals that 

L 

appeared in (13): 

-L z . ~L 

j dz j (ieA(e)sin[g(^-e)] = e-^''^^Res[A(e)e-^''«] + - / A{z)dz 
j dzj rfCA(Oe^''«sin[g(z-C)] = / ^{z)dz 

L L L 

These expressions are derived by expanding the sine into components, integrating by 
parts, retaining only the exponents that increase in the large semicircle, and proceeding to 
integration over the hole large circle as when deriving (18)-(21). Taking into account that 
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(integration goes over the large semicircle), it is easy to obtain the following final formula 
for the absorption spectrum: 



4. Discussion 

The first term in (23) yields a singlet and does not depend on the presence of defects. The 
second term yields a wing in the exciton zone. It is somewhat unexpected that its shape 
does not depend on the relative positions Rr of the defects. The wing amplitude turns out to 
be proportional to the sum of the amplitudes of the defects; this calculation has a "physical" 
level of rigor. For a rigorous calculation , one should write the entire perturbation series in 
A for the solution of (4) and make sure that the infiuence of all terms of this series, except 
the above-considered first-order term, tends to zero in the large semicircle as L — ^ oo. This 
was done and result (23) turned out to be valid only for real q, i.e., for energies falling 
within the excitone zone of the defectless chain. For other energy values, the contribution 
from terms of orders higher than 1 can not be ignored. This behavior should be discussed. 

As shown in@, '^{z) is related to the Green's function of the system and, thus, solution (23) 
could be analytically extended to the range of energies beyond excitone zone of the defectless 
chain. This would be correct if we were dealing with exact Green's function. However, Eq. (4) 
is derived in@ in the continual approximation. Therefore, (23) can yield noticeable deviations 
from the exact spectrum beyond the excitone zone, where the calculation of the previous 
section is invalid. Computer analysis confirms this remark. 

Figure 1 (there is no figure in this version of the paper, but it is easy to imagine) displays 




(23) 



E = n + i6 



the logarithmic absorption spectra of two chains with different numbers, positions and am- 
phtudes of the hyperbohc defects at W — 1, p = 1,5 = 0.02, R — 2. In all cases the number 
of particles was 600, i.e., 2L = 600. The upper part a shows the logarithmic absorption of 
a chain with five hyperbolic defects. Their positions (numbers of the sites at which atomic 
fission becomes infinite, i.e., these atoms are effectively absent) are 200,250,300,350,400. 
The respective amplitudes of the defects are 2, 3, 3, 6, 1. The lower part b shows the same 
for a chain with two defects having the amplitudes - 5 and 5 and positions - 250 and 350. 
For each numeric spectrum, a spectrum calculated by (23) (smoothed curves) is given. No 
fitting (either in amplitude or in shape) was carried out. The energy plotted in units of 
W. It is evidend from Fig.l that for < Q < 1 (the excitone zone), the numerical spectra 
almost completely coincide with the theoretical spectra and the shape of the spectra does 
not depend on the size or positions of the defects. The discrepancy for f2 fti 1 is due to 
the finiteness of the chain because, in this energy range, g^(g^ ~ 0) can be comparable to 
A(z), which is finite in the large semicircle due to the finiteness of the chain. In the region 
beyond the exciton zone, only agreement "in the mean" is possible (it can be improved by 
increasing the imaginary part of the energy S). In this region, the spectrum does depend on 
the positions and sizes of the defects. It may be possible to apply the proposed approach to 
more realistic models of chains with defects. 
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